The Geometry of Crumpled Paper 
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We measure the geometry of a crumpled sheet of paper with laser-aided topography and discuss 
its statistical properties. The curvature of an elasto-plastic fold scales linearly with applied force. 
The curvature distribution follows an exponential form with regions of high curvature localized 
along ridges. The measured ridge length distribution is consistent with a hierarchical model for 
ridge breaking during crumpling. A large fraction of the ridges are observed to terminate without 
bifurcating and the ridge network connectedness is not as complete as anticipated. The self-affinity 
of the surface is characterized by a Hurst exponent of 0.72 ± 0.01 in contrast with previous results. 



A crumpled piece of paper is an interesting and ubiqui- 
tous example of a stress induced morphological transfor- 
mation in thin sheets. As a sheet is deformed, the bend- 
ing energy becomes localized and the resulting crumpled 
shape is often thought of as a network of connected line- 
like ridges. These characteristic topologies are found in a 
variety of objects ranging from biological systems to en- 
gineering applications - for example, polymerized vesicle 
membranes and crumple zones in automobile bodies. In 
all situations, the response of the material can vary from 
purely elastic deformations to elasto-plastic depending 
on the applied stress. Much of the theory for crumpled 
membranes and shells focus on either equilibrium trans- 
formations found in biology 0, or the elastic properties 
of macroscopic sheets |^ Q ■ A number of recent stud- 
ies 0, IE S 0i treat the the geometry of developable 
cones which form the geometrical foundation for struc- 
tures found when elastic sheets are subjected to point 
deformations. However, few studies exist that directly 
explore the crumpled state when elasto-plastic deforma- 
tions occur. 

Early experiments focused on the fractal dimension of 
the ball that circumscribes a flat sheet once crumpled 8j . 
The acoustic emission from sheets that have been crum- 
pled and unfolded display power-law scaling in the dis- 
tribution of click energies ^ J^]. By postulating that 
the ridge lengths are proportional to the energy released 
in acoustic emission, the authors indirectly measured the 
ridge distribution. However, very few experiments which 
directly access the surface leave many fundamental ques- 
tions unanswered. 

In this paper, we report new experiments designed to 
directly measure the geometry of crumpled paper us- 
ing non-invasive, laser-aided topographical reconstruc- 
tion. Large sheets of paper are crumpled to a ball of 
fixed radius and then unfolded to reveal the plastic de- 
formations made by the crumpling process. The surfaces 
are scanned by a laser sheet and the coordinates of the 
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FIG. 1: (a) Schematic diagram of the experimental apparatus, 
(b) Example of a crumpled paper measured with laser-aided 
topography and rendered with shadow lighting. 



surfaces are reconstructed. Thus, we directly access the 
local curvature of the sheet and measure the properties 
of the ridges and their structure. We also examine the 
curvature of a single fold as a function of applied force 
to probe the elastic and plastic deformation which occur 
during crumpling. 

A schematic of the experimental configuration is shown 
in Fig. n^a). Three kinds of paper are used (see Ta- 
ble^) The paper are first crumpled to "spheres" of radius 
i?c ~ 120 mm and then unfolded to a partially flat state. 
The uncrumpling process consists of opening the sphere, 
taking care not to tear the sheet while also making sure 
that we do not flatten the scars that are left on the pa- 
per. The opened sheets are then placed into the scanning 
setup for measurement. A laser with a cylindrical lens 
produces a sheet of light approximately in the vertical di- 
rection and a CCD camera with a resolution of 1024 x 768 
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Mass 


I mm ^ ) 


hp (mm) 


Size (mm x mm) 


1.16 


xlO"'* 


0.156 


560 X 610 


1.89 


xlO"*^ 


0.207 


560 X 610 


2.97 


xlO"'* 


0.360 


430 X 355 



TABLE I: Mass, thickness, and size of paper used in the crum- 
pling experiments. 



pixels is placed at an angle to image the intersection of 
the laser light with the surface. When imaged from the 
side, the bright points where the laser light intersects the 
surface is proportional to its height [see Fig. ^a)]. Each 
column of the image is first scanned for the highest pixel 
value and that element is used for the brightness weighted 
centroid along that particular column. Care is taken to 
ensure that the camera is mounted so that the laser light 
does not "disappear" behind the structure that is being 
scanned. The laser sheet is rotated with a stepper motor 
which allows the entire surface to be scanned. After ap- 
propriate scaling and calibration, the crumpled surface is 
measured to within the thickness of the paper (see [13 
for further details.) An example of a crumpled surface is 
shown in Fig. ^b), and corresponds to a 450 mm x 588 
mm sheet of paper. The resulting "surface" is a matrix 
of elements whose value defines a height. Morphological 
features are extracted using a local quadratic fit to a local 
fitting region 

Before analyzing the entire crumpled surface, it is help- 
ful to examine a single fold and the forces that cause 
them. Therefore, we folded a sheet of paper with thick- 
ness hp = 0.207 mm, width x = 150 mm and length 
y — 588 mm as shown in the inset to Fig. [21 about the 
center of the y-axis (lengthwise). The two edges of the 
folded paper are affixed to the laboratory bench to pre- 
vent slipping, then a known mass (force F) is placed 
across the length of the fold for uniform compression. 
After the mass is placed on the fold the radius of the fold 
i?o is directly measured. After approximately 30 seconds 
of waiting time the mass is removed and the paper is un- 
folded and placed in the laser scanning apparatus. The 
inset to Fig. [3 shows schematically the resulting form of 
the unfolded sheet which is plastically deformed. A least 
squares fit to a circle centered at the apex of the fold 
yields the radius of curvature, Rf. It should be noted 
that the radii are measured on the outer surface and 
therefore include the thickness of the paper used. These 
measurements imply a lower bound on each measured 
radius, Roj > hp where hp is the paper thickness. 

In Fig.[21the initial and final curvatures C, {R^^,RJ^) 
are plotted versus the force per unit length F' — F/y. 
We observe that as expected, the curvature of the fold 
increases with an increasing applied force. It is also in- 
teresting to note that although there is always plastic 
deformation present, the data is linear over two orders of 
magnitude in applied force. Furthermore, the timescales 
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FIG. 2: The curvature C of a folded sheet plotted versus F' , 
the force per unit length. The shapes corresponding to the 
folded (v) and unfolded (o) curvatures are shown in the inset. 
These points give a calibration for the forces that cause the 
deformation in crumpled paper. 



hp = 0.156 mm - 
hp = 0.207 mm- 
hp = 0.360 mm - 




FIG. 3: The probability distribution of the magnitude of the 
curvature P(c) averaged over regions of 25 mm^ for each thick- 
ness of paper. The increase of the curvature for thicker sheets 
indicates that the uncrumpling process does not removed the 
quenched curvature of the ridges for thicker sheets. 



for the application of the applied force in the hand crum- 
pling and the single folding are equivalent allowing for 
direct comparisons of the deformations. 

With this understanding of the forces involved to pro- 
duce a single fold, we now turn to the analysis of a fully 
crumpled sheet. One of the parameters used to assign a 
label to the point of interest is the cross-sectional curva- 
ture c, defined as a plane that intersects with the plane 
of the slope normal, and perpendicular aspect direction. 
The distribution of the curvature magnitudes averaged 
in boxes of size 5 mm^, c, are shown in Fig. |31for all hp. 
Note that as hp increases, the high curvature tails also 
increase, indicating that for thicker sheets, plastic defor- 
mation is well preserved after uncrumpling. This result is 
expected by simply considering that the energy required 
to unfold a ridge must be proportional to hp. 

A proposed model of the distribution of energy in 
a crumpled sheet states that the process of crumpling 
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FIG. 4: The distribution of measured ridge lengths P{£) (a) 
hp = 0.156 mm, (b) hp — 0.207 mm, (c) hp — 0.360 mm 
display. The dashed line is a fit to Eq. Q and shows good 
agreement and captures the existence of a peak, (d) Width of 
the log-normal distribution a versus the paper thickness hp, 
the dashed line is a linear least squares fit. 

forces the majority of the energy imparted to the sheet 
will be localized in a network of line- like ridges 0, Q, 
lisl flil . The structure of the ridge network is said to be 
that of a simple hierarchical structure. The structure is 
composed of large ridges that are randomly and unevenly 
"broken" into smaller ridges that produces a set of frag- 
mented lengths. Therefore, a ridge of length £ should be 
equivalent to the original ridge length times a set of ran- 
dom variables that correspond to the fraction of ridge left 
after each breaking. The logarithm of the ridge lengths, 
becomes a sum of random variables, i.e. a random walk. 
Therefore, the probability density function of the ridge 
length has the form, upon a change of variables back to 
lengths [l^ . 

dl a'^i 

A log-normal distribution with a pre- factor. To test 



this prediction, we measured the lengths of the ridges 
found in the hand crumpled sheets by identifying the 
points that lie along individual ridges via a nearest neigh- 
bor criteria. We then perform a spline fit to each individ- 
ual ridge and calculate the contour length of each spline. 
In Fig. 0] the distribution of ridge lengths are plotted for 
each hp. The dashed line is a best fit to Eq.^ To within 
experimental accuracy, the data is well described by this 
distribution function. 

In theoretical treatments of crumpled surfaces, either 
equilibrium P, 0|, or non-equihbrium 0, 0, 0], the 
structure of the surface is comprised of randomly ori- 
ented ridges that always intersect at vertices to form a 
highly connected network. However, we observe that the 
connectedness of the network is not complete as antici- 
pated. We also measure the number of neighbors that the 
"end" of a ridge will have, and find that many ridges do 
not intersect with other ridges. In Fig. [Sfa), we plot the 
histogram of the number of nearest neighbors iV„, for 20 
individual hp = 0.207 mm sheets. The existence of ridges 
with zero nearest neighbors could have two explanations. 
Elastic deformations do not leave scars of the crumpling 
process. Therefore, the ridges that existed during crum- 
pling do not appear when the sheet is unfolded. Second, 
the paper itself may be able to absorb the force that cre- 
ates the ridge. The ridge may essentially dissipate at it's 
edges due to the thickness of the paper itself. A question 
may arise if the sheet is crumpled "hard" enough. To 
answer this question, we note that squeezing the paper 
into a smaller Rc than used here often results in tearing 
which introduces new complications. 

We also measure the angles that the ridges create with 
their neighbors. In Fig.|3Jb) the histogram of the angles 
6, that are created by four ridges is plotted {i.e. the data 
in the iV„=3 bin which corresponds to single d-cones), 
when they meet at a vertex. The values of 9 appear 
broadly distributed in the range > > 180 deg. indi- 
cating a somewhat random ordering. However, a number 
of predominant peaks at < 20 deg. « 60 deg. and 9 w 
110 deg. are apparent. The existence of these peaks in- 
dicates that indeed the d-cone geometry dominates the 
intersection of ridges when Nn — 3 consistent with pre- 
vious results 0,1131. However, the precise value of the 
predominant opening angle is still debatable. 

Finally, to draw connections with previous experimen- 
tal and theoretical studies of crumpled surfaces [ij, we 
also performed surface scaling measurements. A self- 
affine surface C,{x,y), is defined by the scaling of that 
surface by the transformation 

X ^ Ax, y ^ Ay, C ^ A^C, (2) 

where A is an arbitrary scaling factor and H is known 
as the Hurst exponent. We measure H directly for each 
dimension using the following method. Each line scan 
is sectioned into non-overlapping bins of width w that 
are increased in size until w < L/2, where L is the 
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to be consistent with universal scaling laws, many sub- 
tleties are also revealed. The existence of a linear re- 
lationship between the applied stress and the resultant 
curvature of the plastic deformations implies that the sur- 
face curvatures can be mapped to the force distribution of 
crumpled surfaces. Due to the linearity of the response, 
a simpler theory may be possible to explain our observa- 
tions. Additionally, we observe a considerable fraction of 
ridges end without bifurcating. Some of the differences 
may arise because of the fact that virtual surfaces which 
can pass through each were considered in previous analy- 
sis. Including interactions between surfaces as the paper 
is squeezed presents new and interesting challenges that 
needs to be addressed by further theoretical work. 

We thank T. Witten, R. Mukhopadhyay, L. Mahade- 
van, and A. Boudaoud for stimulating discussions. 



FIG. 5: (a) Histogram of the number of nearest neighbor 
ridges at a vertex Nn . The high frequency of zero occurrences 
indicates that many ridges dissipate into the paper, (b) His- 
togram of angle between ridges 9 for A'^„ = 3. (c) Hurst plot 
of the one-dimensional line crossections of the crumpled pa- 
per. The least square fit gives the Hurst exponent H. We 
observe that for small w, there exists a scaling usually asso- 
ciated with 1// noise and that for larger intervals the scal- 
ing follows H ~ 0.71. (d) Hurst plot of the two-dimensional 
crumpled gives H ^ 0.72. {hp = 0.207 mm) 



length of the scan. To ensure that we may consider each 
scan as statistically independent, we skip every ten scans. 
For each bin, the maximum difference between elevations 
within that bin A, is recorded. This process is continued 
for each line and the average of each interval is taken. For 
two dimensions the same process is used. The width w 
now represents the length of a box where the maximum 
difference is measured. The relationship. 



(3) 



directly gives the Hurst exponent. The one-dimensional 
Hurst plot [see Fig. |SJc)] demonstrates two distinct scal- 
ing regions. For small w, H 0.95 indicating that 
for very short intervals the data seems to be consis- 
tent with 1// noise. For larger intervals, the scaling 
follows H ^ 0.71. The two-dimensional analysis [see 
Fig. |SJd)] also demonstrates Hurst scaling. Utilizing 
one-dimensional scanning profilometry of crumpled pa- 
per Plouraboue and Roux 21] have reported H = 0.88. 
In a theoretical investigation, based on a lattice model 
of crumpled paper, Tzschichholz, et al. reported 
H = 1.0. Both values are in contrast with our results. 

In summary, we have examined the structure of a 
crumpled sheet of paper using laser aided surface topog- 
raphy that allows for the complete reconstruction of the 
crumpled surface to unprecedented precision. While cer- 
tain features such as ridge length distributions are found 
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